Electronic compressibility of layer polarized bilayer graphene 
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We report on a capacitance study of dual gated bilayer graphene. The measured capacitance 
allows us to probe the electronic compressibility as a function of carrier density, temperature, and 
applied perpendicular electrical displacement D. As a band gap is induced with increasing D, the 
compressibility minimum at charge neutrality becomes deeper but remains finite, suggesting the 
presence of localized states within the energy gap. Temperature dependent capacitance measure- 
ments show that compressibility is sensitive to the intrinsic band gap. For large displacements, 
an additional peak appears in the compressibility as a function of density, corresponding to the 
presence of a 1-dimensional van Hove singularity (vHs) at the band edge arising from the quartic 
bilayer graphene band structure. For D > 0, the additional peak is observed only for electrons, 
while D < the peak appears only for holes. This asymmetry that can be understood in terms of 
the finite interlayer separation and may be useful as a direct probe of the layer polarization. 



The unique band structures of monolayer (MLG) and 
bilayer graphen^^^ (BLG) offer unprecedented tunabil- 
ity in a high quality two dimensional electron system 
(2DES). Within the independent electron approximation, 
both MLG and BLG are gapless, chiral systems. The 
gapless spectra are related to the pseudospin degeneracy, 
which is tied to the symmetry between the two sublat- 
tices constituting the honeycomb. Whereas in the mono- 
layer a gap can be opened only by a potential modula- 
tion on the spatial scale of the lattice constant^, in BLG 
the relevant sublattices are located on different layers, 
allowing a gap to be induced by a modulation of the in- 
terlayer imbalance via the application of an electric field 
perpendicular to the BLG planes^^. Although the field- 
effect tunable gap in BLG has been observed optically 
^^^^ transport measurements show hopping conductivity 
at low temperatures.'^^"^^' 

In a parallel plate capacitor made up of imperfect con- 
ductors, adding charge n to the plates costs the sum of 
the classical electrostatic energy, the kinetic energy due 
to the resulting change in the chemical potential /i, and 
the potential energy of Coulomb interactions between the 
charge carriers. The measured differential capacitance, 
C = Sn/SV^ in such a system reflects this finite electronic 
compressibility by manifesting a lowered effective capaci- 
tance, = Co~^ + (e^Ai^) , where Cq is the geometric 
capacitance, A is the area of the device and v = dn/djii^ 
the electronic compressibility, which corresponds to the 
density of states in the noninteracting, zero temperature 
limit. In low dimensional systems the contribution of 
the compressibility to the capacitance — termed "quan- 
tum capacitance"^^ — can be small even when the con- 
ductivity remains large, providing a powerful tool in the 
study of both one^^ and two-dimensional electronic sys- 
tems. Capacitance measurements are particularly power- 
ful in the study of disordered systems, as they are able to 
detect localized states whose contribution to transport is 



suppressed. Capacitance measurements are , as a result, 
crucial in understanding biased bilayer graphene, a sys- 
tem in which localization is known to play a role. More- 
over, the small but finite interlayer separation {d ~ 3.4 A) 
allows the layer polarization to be probed through elec- 
trostatic measurements, as described in a companion pa- 
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per"^' . 

To produce dual gated graphene devices with high ge- 
ometric capacitance, we utilize single crystal hexagonal 
boron nitride (h-BN) flakes^^ as the top gate dielectric 
fabricated by the process described in Ref. (191 Briefly, 
both graphene and single crystal h-BN, an insulating iso- 
morph of graphite, are exfoliated onto n-Si/Si02 wafers 
(Fig. [ija-b)). A thin (5-7 nm) h-BN flake is trans- 
ferred on top of the graphene usinga wet etch process 
and micromechanical manipulation^, followed by elec- 
tron beam lithography to form contact electrodes and 
a local top gate (Fig. [l]^). The heavily doped silicon 
substrate, coated with 285 nm oxide, serves as the bot- 
tom gate. The double gated geometry allows independent 
control of the electronic density, n, and the displacement, 

S = H ((Fb - - g| {Vt - F^)) through the top 
and bottom gate voltages, Vt and Vb- Here, Eb and 
dB are the dielectric constant and thickness of the back 
gate dielectric layer, C^(^^^ is the geometric capacitance 
of the top (bottom) gate, and V^^^^^ is the voltage offset 
required to obtain minimal density and displacement in 
the dual-gated region. We find that h-BN is an excel- 
lent gate dielectric, with e ~ 3-4 and breakdown fields 
comparable (~.8 V/nm) to Si02 thin films. In addition, 
we observe minimal degradation of graphene samples, 
with no additional doping contributed by the presence 
of the top gate and typical post h-BN transfer mobilities 
of /i ~5, 000- 10, 000 cm^/V sec for graphene monolayers 
and /i ~2, 000-3, 000 cm^/V sec for bilayers. 

Low temperature capacitance measurements were per- 
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FIG. 1: (a-c) Optical microscope images, (a) a thin hBN sin- 
gle crystal is transferred onto a (b) mono- or bilayer graphene 
flake and (c) contacts and gate electrodes deposited by elec- 
tron beam lithography. Scale bars are 10 /xm. (d) Cross sec- 
tion schematic of resulting dual gated device, (e) Schematic 
circuit diagram of the capacitance bridge. A reference capac- 
itor (Johanson Technology R14S) is mounted on the probe, 
and a reference voltage is chosen to balance the capacitance 
bridge. A small AC excitation signal is added to the DC gate 
bias through a transformer (Triad Magnetics SP67). 



formed using a capacitance bridge circuit (Fig. 1(e)) 
with a cold reference capacitor. All wires were shielded, 
and the sample package was encased in a Faraday cage 
to further reduce parasitic capacitances, which represent 
an additive constant to the measured value of Ct- A 
ceramic multilayer capacitor with minimal temperature 
dependence was chosen for the reference capacitor {Cr) 
and connected near the sample at low temperature. The 
noise level of the bridge was ~25 e/Vi^z, allowing sub- 
femtofarad resolution with averaging times of less than 
30 seconds for our typical top gate AC excitation volt- 
age 5Vt ~ 15-50 mV. Measurements were performed off- 
balance. For the variable temperature capacitance mea- 
surements, we applied 5Vt =50 mV AC excitation volt- 
age on top of the DC gate bias and measured the cur- 
rent through the graphene device directly. Although this 
method results in poorer signal to noise than the bridge 
measurement, it eliminates calibration errors stemming 
from the small temperature dependence of the reference 
capacitor. The ability to measure the quantum contribu- 
tions to the capacitance relies on the use of a thin top gate 
dielectric layer. The 7 nm thick hBN dielectric results in 
values of that are comparable to the quantum capac- 
itance Cq, so that variations in the measured Ct with 
changing electronic compressibility are easily detectable. 
We check the efficacy of this measurement scheme us- 



FIG. 2: Measured conductance and capacitance of monolayer 
graphene at 2 K and B = (a) and B = 9 T (b). In (a), ca- 
pacitance and conductance resemble each other closely due to 
the fact that both contain a spurious contribution adding "in 
parallel"; G = (1/i? + l/i?c)"^ and Ct = {I/Cq + l/C^)'^ 
where R is the device resistance, Cq is the compressibility, 
or "quantum capacitance," C^ is the geometric capacitance, 
and Rc is the contact resistance of our two-terminal devices. 
In (b), capacitance reflects the electron 



ing monolayer graphene capacitors. MLG is expected to 
display monotonically increasing as a function of ab- 
solute density \n\. In the presence of a strong magnetic 
field, V is further modulated due to the formation of Lan- 
dau levels. Fig. [2] shows the measured capacitance C and 
conductance G of MLG as a function of top gate voltage 
Vt at both zero and finite magnetic field. The lowered 
compressibility stemming from the linear spectrum of of 
MLG can be inferred from a depression in the capacitance 
at zero density, while at high B the formation of the zero 
Landau level leads to a peak at charge neutralitjEU. The 
high magnetic field capacitance traces show compressibil- 
ity oscillations due to the formation of higher LLs, while 
conductance shows the electron-hole asymmetry that is 
the signature of edge state transport in graphene hetero- 
junction^^^H^, which form due to the partial coverage of 
the graphene channel by the top gate. The peak inver- 
sion at zero density and electron-hole symmetry in high 
magnetic field together confirm that the measured signal 
is indeed capacitive, and is sensitive to density of states 
rather than lateral resistance effects. 

We now turn our attention to bilayer graphene (BLG) 
samples. Fig. |3] shows the capacitance of a BLG sample 
at 1.5 K measured with the cold bridge. Tuning external 
gates adjusts both n and D. For small values of I) ^ 0, 
the measured capacitance exhibits a minimum at n = 
as expected for ungapped bilayer graphene, which has a 
hyperbolic band structurd ^^ l ^^ l As increases, the n = 
minimum gets deeper, corresponding to the formation 
of a gap in the energy spectrum. The n=0 minimum 
does not go to zero at high values of \D\^ and in fact 
the capacitance modulation is only 10% with respect to 
the D = {) value. In addition, a distinct local maximum 
develops next to the minimum. As we argue below, the 
presence of both the dip and local peak in v at high \D\ 
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FIG. 3: Capacitance at B=0 and 1.5 K as a function of Vt and Yb- Colored traces in (b) are taken at 30 volt intervals in Vs, 
corresponding to the colored lines in (a), (c) Traces at constant D, extracted from the data set shown in (a). Data is plotted as 
a function of C^^AVr + C^AVb, which would correspond to the density were the bilayer perfectly 2 dimensional and perfectly 
compressible. Curves in (c) are offset for clarity by 2fF per V/nm in D. 



can be understood, at least qualitatively, from the band 
structure of gapped BLG, once the effects of disordei!^ 
and the interlayer separation^^ are taken into account. 

Within the nearest-neighbor tight binding approxima- 
tion, the energy spectrum of pristine, Bernal stacked 
bilayer graphene with finite interlayer asymmetry A is 
gapped and has a "Mexican hat" structure^^. Even 
in the presence of disorder, the absence of a positive 
quadratic term in the energy spectrum turns the problem 
of gapped, disordered bilayer graphene into one of a heav- 
ily doped semiconductor with quartic energy band^. 
We thus expect vestiges of a ex 1/ vHs-like feature 
to be present even in our low mobility bilayer samples, 
manifesti ng as a nonmonotonic-in-density feature at the 
band edg^^^. In addition to smearing the band edge 
vHs, disorder has a dramatic effect on compressibility at 
charge neutrality in the presence of a large gap. In con- 
trast to clean semiconductors, in which the depleted sys- 
tem is incompressible, our measured capacitance remains 
finite and large even for large a fact we attribute to 
tails in the density of states representing localized in- 
tergap stated. This explains the discrepancy between 
energy scales that govern transportP^HI^ and the gap en- 
ergies observed opticall}CH^. It also suggests that recently 
predicted topological edge conduction^^ is not the domi- 
nant reason for incomplete "turn-off" in BLG devices of 
typical quality. 

Quantitative analysis of the capacitance data requires 
extracting the compressibility u from the measured sig- 
nal. For a perfectly two dimensional electron system 
{Cbl ^ oo), the measured top gate capacitance is 

where Cpara are all parasitic capacitances between gate 
and contact electrodes and terms O (^^f^ ~ -034 have 



been neglected. Extracting u thus requires subtract- 
ing both parallel (Cpara) and series (C^) capacitances, 
and dividing by A =31/im^ (determined by optical mi- 
croscopy). We determined C^/C^ from the ratio of the 
back and top gate capacitances, measured by tracking the 
charge neutrality point in the Vt-Vb plane (the dark 
diagonal belt in Fig. |3ja)). For the BLG device pre- 
sented in this paper, we measured and Ct/Cb =29.5±.4, 
where Cs = 115 aF//im^ is the geometric capacitance of 
the bottom gate. As disordered BLG devices cannot be 
turned off completely, Cpara cannot be measured in situ 
as is common practice in depletable semiconductors 
and semiconducting carbon nanotubes^^. Instead, we de- 
termine Cp=16±l fF by removing the graphene through 
a short oxygen plasma etch and measuring the remain- 
ing capacitance between the metal contacts and the top 
gate. Cpara coustitutcs about 10% of the total capaci- 
tance signal. Due to the subsequent subtraction of the 
(series) geometric capacitance C^, the error in determi- 
nation of Cpara IS Icast important when the capacitance 
differs considerably from the geometric value. This is the 
regime in which we perform a quantitative analysis of the 
compressibility of gapped BLG. 

Near overall charge neutrality at \D\ ^ our samples 
show a hopping conductivity similar to that observed in 
Ref. [To] from 1 K to ~ 150 K. Capacitance instead shows 
no significant temperature dependence up to 50 K, there- 
after slowly rising as might be expected for a spectrum 
with a gap in 50-100 meV range (Fig. [4]). This is con- 
sistent with the presence of disorder-induced tails in the 
density of states throughout the band gap: whereas tem- 
perature dependent transport is dominated by the hop- 
ping between these localized states, temperature depen- 
dent capacitance is dominated by thermal population of 
the much larger density of states near the band edge. 

The most interesting, and unexpected feature, of the 
experimental data are the local maxima observed at the 
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FIG. 4: Bottom panel: Temperature dependence of the in- 
verse compressibility (AeV)"^ for Vb = 80, and -80 V. 
A single value of Cpara is chosen for all gate voltages at a 
given temperature, but a different Cpara is chosen for each 
temperature so that the curves match at high density. We 
attribute the variation in Cpara to thermal expansion of the 
bonding wires, the capacitance of which constitute the major- 
ity of Cpara- Top panel: comparison between the temperature 
dependence of the minimal compressibility and minimal con- 
ductivity at Vs=-60 V. The dashed line is a guide for the eye. 



band edge are associated with the ID vHs inherent in 
the BLG band structure EllMI. Interestingly, this fea- 
ture is only present on one edge of the band, appear- 
ing on the electron side for D > and the hole side 
for D < 0. This inversion symmetry in the variables 
(D^n) was observed in all devices measured, including 
those fabricated using a resist free shadow mask met- 
allization as well as seeded atomic layer deposition of 
Hf02^^. Understanding this asymmetry requires taking 
into account the three dimensional structure of BLGP^, 
which consists of two strongly coupled but spatially dis- 
tinct layers of carbon atoms. The charge distribution on 
a BLG flake is sharply concentrated on the two layers, 
n{z) :^ ni6{z) + n2S{z + (i), so that the system can be 
modeled as a four plate capacitor. Solution of this elec- 
trostatic problem leads to the modified relation for the 
measured capacitance in which Ct depends on both the 
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FIG. 5: Left Panel: Subtracted capacitance, ACt = Ct{D) — 
Ct{D = 0), as a function of approximate density for different 
applied displacements. Curves are offset by IfF per V/nm in 
D. Right Panel: Calculated top gate capacitance for disor- 
dered bilayer graphene, following Ref. 17 The colors cor- 
respond to the displacements in the left panel; the disorder 
parameter is 7 = 4 for all curves. Curves are offset by 3 fF 
per V/nm in D. 



interlayer capacitance, V2i^ and intralayer capacitances, 
vii and 1^22, where Uij = dui/dvj with n and v rep- 
resent layer indexed density and potential, and the in- 
dices i,j=l,2 denote the top or bottom layer. Crucially, 
while V21 is symmetric with respect to layer interchange, 
1^ is obviously asymmetric in the presence of interlayer 
asymmetry. Penetration field^^ measurements of bilayer 
graphene depend only on layer-symmetric quantities, and 
thus probe fundamentally different physical quantities. 

As elaborated in Ref. [17| the 1 / y/E divergence associ- 
ated with the ID vHs at the band edge principally on the 
low energy layer within the BLG flake. The vHs man- 
ifests more strongly in the measured capacitance when 
vHs- hosting layer is closest to the top gate; conversely, 
the gate sees a far-layer vHs only through the screened 
field penetrating the near layer. Counterintuitively, dis- 
order enhances this effec t not only by smearing the to- 
tal density of state^^^^^ but by populating the normally 
depleted non-vHs layer, thereby enhancing its ability to 
screen. While an ideal experimental geometry would per- 
mit the simultaneous measurement of capacitance from 
two sides of the BLG flake this is effectively accomplished 
in our single local gate geometry by reversing the sign of 
thus reversing the order of the vHs bearing and non- 
bearing layers. 

In order to better compare our experimental data with 
theory based on a parabolic two band model, it is conve- 
nient to subtract a background taken at I) = 0. Because 
the high energy behavior depends only weakly on the 
displacement, this has the effect of isolating the low en- 
ergy part of the measured capacitance, and in addition 
removes the effects of electron-hole asymmetric elements 
of the band structur^^. The results of this subtraction 
resemble theoretical calculations which take into account 
both the interlayer separation as well as weak, short- 
range disorder (Fig. [5|. In particular, the asymmetric 
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appearance of the van Hove singularity can be under- 
stood as the effect of disorder enhanced interlayer screen- 
ing. However, quantitative understanding of the role of 
disorder in bilayer graphene will require experiments that 
independently control the disorder, as well as a more so- 
phisticated theory taking into account a wider variety 
of effects including long range scattering and electron- 
electron interactions. 

In conclusion, we study the broken symmetry state of 
bilayer graphene induced by an electrical displacement 
field applied perpendicular to the bilayer. We observed 
the formation of the displacement induced gap and its ac- 
companying ID vHs, and estimate the gap size from tem- 
perature dependent measurements. We also discussed 
how different measurement geometries make capacitance 
a probe of layer-pseudospin polarization in the bilayer. 
In BLG, layer symmetry breaking in BLG can also be 



expected to occur spontaneously, as a result of electronic 
interactions at both finite^^ 33 36 zero^^ magnetic 
fields. Similar measurement performed on high-quality 
bilayers manifesting these effects will allow spontaneous 
layer polarization to be probed directly through compar- 
ison of measurements in various capacitance geometries, 
as discussed at length in a companion papei^-^. 
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